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The variances and covariances associated to the nuclear matrix elements (NME) of neutrino- 
less double beta decay (OV/3/3) are estimated within the quasiparticle random phase approximation 

■ (QRPA). It is shown that correlated NME uncertainties play an important role in the comparison 

of 0vf3f3 decay rates for different nuclei, and that they are degenerate with the uncertainty in the 

, reconstructed Majorana neutrino mass. 
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I. INTRODUCTION 

The search for the neutrinoless mode of double beta decay (0vf3f3), 

Oh" {Z, A) -> (Z + 2, A) + 2e- , (1) 

i ' 

is being vigorously pursued by several experiments using different (Z, A) nuclei, in order to unravel the Dirac or 
Majorana nature of neutrinos and their absolute mass scale 1]. In a given candidate nucleus i = (Z, A), light 
i— > Majorana neutrinos can induce 0v/3f3 decay with half-life T, given by 
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■ where Gi is a calculable phase-space factor, M[ is the 0z//?/3 nuclear matrix element (NME), and mpp is the "effective 
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Majorana neutrino mass," 
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where and U e k are the neutrino masses and the v e mixing matrix elements, respectively, in standard notation [2j. 
The NME includes both Fermi (F) and Gamow- Teller (GT) transitions, plus a small tensor (T) contribution 0], 

Si «M^r( M - T+M --f)- w 

In the above expression, qa is the effective axial coupling in nuclear matter, not necessarily equal to its "bare" free- 
nucleon value qa — 1-25. With the conventional prefactor oc g\ in Eq. ((4]), the phase space Gj becomes ^-independent. 
In general, all parametric uncertainties (which may be quite large) are embedded in \M[\ 

It is widely recognized that a convincing case for CV/3/3 decay must involve independent signals in three or more 
nuclei [J [5[. For instance, if the theoretical NME uncertainties could be roughly expressed in terms of a single 
nuisance parameter p, then one would need two independent half- life data T% and T2, and two relations as Eq. I[2]). to 
fix both p and mpp (up to degeneracies). A third datum T 3 would overconstrain the system of equations, providing a 
cross-check of the results 0, 0, 0- A negative check might signal possible new Qv{30 physics beyond light Majorana 
neutrinos (barring experimental or theoretical mistakes). Any new 0v/3f3 mechanism(s) would then involve at least one 
more unknown, and thus it might require one or more data (T4, T5, . . .) for further cross-checks (9I. Iiol. . Statistical 
assessments of the various options demand realistic estimates of experimental and theoretical uncertainties, and the 
analysis of possible degeneracies which, as we shall see, may play a relevant role. 

Recently, there has been significant progress towards the reduction (and a better evaluation) of 0v(3f3 theoretical 
errors. Within the quasiparticle random phase approximation (QRPA) Q, these uncertainties can be largely kept 
under control by systematically fixing, in each nucleus, the particle-particle strength parameter g pp via two-neutrino 
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double beta (2vj3j3) decay rates. In this way, the dispersion of NME values obtained by varying several QRPA 
ingredients has been significantly reduced (see [J] and references therein). 

However, the estimated NME variances do not exhaust the information needed to compare 0v(3(3 limits (or signals) 
in different nuclei: the NME covariances are important as well. Nonzero NME covariances have been implicitly 
recognized in a few works, e.g., by studying the dispersion of NME ratios [81], and by observing that such a disper- 
sion may be smaller than for individual NME [io| . To our knowledge, these observations — implying positive NME 
correlations — have not yet been sharpened from a statistical viewpoint, despite their relevant consequences for the 
comparison of Qvf3(3 signals. In a nutshell, the main points can be illustrated as follows. If a finite half- life Tj is 
measured in a nucleus i, the half- life expected in another nucleus j is 

_ rp Gi \ M i\ 2 (rs 

3 ~ 1 Gj |m;p ' ( j 

within (large) NME uncertainties. From the experimental viewpoint, the "most favorable case" would entail the 
shortest decay timescale Tj, namely, the smallest \M' i \ and the largest |Mj|. However, if the two NME uncertainties 
were positively correlated (e.g., via a common normalization factor), opposite changes of \M[\ and \M'A would be 
unlikely, thus preventing the occurrence of the "experimentally favorable" case. Moreover, a common shift of the 
NME for all nuclei could always be compensated by an inverse shift in rnpp via Eq. ([2]), leading to a degeneracy 
between (correlated) theoretical errors and the Majorana neutrino mass. 

The purpose of this paper is to explore and discuss these issues in detail. In Sec II we set our notation and 
conventions. In Sec. Ill we present our evaluation of the covariance matrix for the NME in a set of nuclei. In Sec. IV 
we apply our formalism to relevant cases in the 0v/3f3 phenomenology. In Sec. V we summarize our work and discuss 
future perspectives. An Appendix collects additional details about different theoretical evaluations of Gj and \M-\. 



II. NOTATION AND CONVENTIONS 

In the spirit of Refs. [3 12111, we shall use logarithms of the main Qvj3f3 quantities in appropriate units, namely: 

n = log 10 (Tj/y) , (6) 
- 7 < = log^GiAy^eV- 2 )] , (7) 

m = io gl0 K| , (8) 

M = logi ("WeV) , (9) 



so that Eq. © is linearized as 



Central values and errors will be denoted as 



t 1=1i - 2rn - 2 M . (10) 



n = t? ± Sl , (11) 

m = v° ± <r< , (12) 

H = fj° ±6 , (13) 

the 7i having virtually no uncertainties (see however the Appendix for remarks) . Experimental measurements of the 
Tj's are thus translated into linear constraints on the unobservable quantity /i, once the nuclear matrix elements r\i 
and their covariances are given. 

Linearization through logarithms is appropriate to deal with relatively large NME errors. For instance, a typical 
"factor of two" uncertainty, \M[\ = \M'®\ x (l^o"), entails at least two drawbacks: (i) asymmetric errors are difficult 
to manage with usual statistical tools (such as least-squares methods); (ii) the unphysical region \M[\ < is hit at 
twice the lower error. Both drawbacks are avoided by expressing the same "factor of two" uncertainty as rji = ^±0.30. 

Concer ning the quantities Tj = log 10 (T,/y), at present there is only one claim for a positive Ov0B result by Klapdor 



et al. [15|,ll6| as part of the Heidelberg-Moscow Collaboration: Tj/y = 2.23±g;^ x 10 25 at la il6j]. We translate this 
claimed range as 

Tj = 25.355 ± 0.072 (la, i = 76 Ge) , (14) 



where we have slightly displaced the experimental central value so as to reproduce the la extrema, by construction, 

l-n-0.44 
-ti-0.31 



with symmetric errors |17|. For n > 1, the asymmetric na ranges Tj/y = 2.23^"'o'gi x 10 25 correspond to the 
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TABLE I: For each nucleus i, we report the phase space factor 74, the central value of the nuclear matrix error r\i, and the 
error <7j, together with the (symmetric) error correlation matrix pij, according to the QRPA estimates in this work. See the 
text for definitions. 
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7i 


V° 
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82 Se 


96 Zr 


correlation matrix pij 

100 Mo U6 Cd 
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100 Mo 


24.660 


0.503 


0.162 


0.973 


0.957 


0.491 


1 










116 Cd 


24.622 


0.404 


0.150 


0.961 


0.961 


0.474 


0.965 


1 








128 Te 


26.073 


0.534 


0.154 


0.947 


0.968 


0.515 


0.916 


0.930 


1 






130 Te 


24.674 


0.498 


0.158 


0.899 


0.927 


0.575 


0.862 


0.870 


0.964 


1 




136 Xe 


24.644 


0.254 


0.187 


0.805 


0.846 


0.663 


0.747 


0.773 


0.898 


0.916 


1 



symmetric ranges r, = 25.355 ± n ■ 0.072 within an acceptable accuracy of 10%, i.e., within about 0.2cr (0.3<r) at the 
level of 2cT (3a) ranges. 

The above arguments, as well as the advantages of using linear relations [Eq. (fTTJ)) ] and the associated simple 
statistics (linear propagation of errors, \ 2 methods), lead us to assume approximately gaussian errors on log Z{, rather 
than on Zi (where z% = 2i, |A£^j), for the purposes of this work. In the future, a better knowledge of the probability 
distributions of the z^s might warrant a different approach, possibly based on more refined statistical tools applicable 
to generic random variables (maximum likelihood methods, MonteCarlo simulations). However, our main results do 
not crucially depend on these subtle aspects. 

III. NME UNCERTAINTIES AND THEIR CORRELATIONS 

In this Section we discuss estimates for the nuclear matrix elements r\i , in terms of central values 77? , errors Oi and 
correlations py, for a set of eight 0i//?/3 candidate nuclei: i = 76 Ge, 82 Se, 96 Zr, 100 Mo, 116 Cd, 128 Te, 130 Te, and 136 Xe. 
We remind that the associated covariance matrix is cov(?yi, r/j) = pijdi/jj, whose diagonal elements coincide with the 
variances of. 

A. Numerical evaluation of QRPA uncertainties 

Our estimates are based on a large set of QRPA calculations 0, [l8[ which include 2x2x3x2 = 24 variants in 
the input ingredients, namely: (i) two values for the axial coupling: gA — 1-25 (bare) and gA — 1-00 (quenched); 
(ii) two approaches to short-range correlations (s.r.c): the so-called Jastrow-type s.r.c, and the unitary correlation 
operator method (UCOM); (m) three sizes for the model basis: small, intermediate and large; (iv) two many-body 
models: QRPA and its renormalized version (RQRPA). All the 24 variants are supplemented by errors induced by 
g pp uncertainties (within the_experimental 2vf3(3 constraints). Concerning NME error estimates, we adopt the same 
conservative approach as in 18], and define the lo range rfl ± at as the one embracing the minimum and maximum 
calculated value of rji for each nucleus i. These ±lcr errors are more generous than their formal statistical definition 
(which would embrace only ~ 68% QRPA variants, i.e., ~ 16 out of 24). Finally, we calculate the correlation index 
between joint (r]i,r]j) values taken from the same QRPA sample. In all cases, we also include <7 pp -induced variations. 

Our final results for 7]i, <Ji and pij are reported in Table I (together with the values of the phase space factors 
7i, for completeness). Figure 1 shows the same results in graphical form, for each couple of different nuclei, in the 
plane charted by the coordinates (rn,r)j). In each panel we show the "lc error ellipse," centered at (rj°,r]j) and with 
correlation pij\ its projection onto a coordinate axis coincide with the ±loi range defined previously. Also shown in 
each panel is the set of QRPA calculations used, supplemented by the horizontal and vertical error bars induced by 
g pp uncertainties (for a total of 24 "crosses" in each plane) . 

In Fig. 1, the strong, positive correlation among theoretical estimates emerges at a glance. The QRPA calculations 
are mostly scattered along a primary direction (the ellipse major axis) with positive slope, essentially as a result of 
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variations in the s.r.c. model (either Jastrow, blue, or UCOM, red) and, secondarily, to variations in gA- There is also 
some dispersion in the orthogonal direction (ellipse minor axis), which is mainly due to g pp variations. In general, the 
overall scatter of QRPA is very well captured by the ellipses, with the possible exception of those involving j — 96 Zr, 
which are somewhat under-sampled at low rjj. For this nucleus, the g pp parameter turns out to be extremely close to 
the so-called QRPA collapse point, the rjj estimates becoming less reliable and more erratic as collapse is approached — 
leading to large and asymmetric error bars. For other nuclei, g pp is far from the collapse point and the results are 
more stable (with smaller and more symmetric g pp errors), as compared to 96 Zr. In conclusion, the correlations p y - 
reported in Table I appear adequate to characterize the scatter of QRPA variants, with the only possible exception 
of 96 Zr, whose estimates must be taken with a grain of salt. 

We remark that the above estimates, performed within the QRPA, include only known and controllable sources of 
uncertainties. Some of them are peculiar of QRPA (e.g., g pp ), while others are common to any nuclear model (e.g., 
gA and the s.r.c). It is not excluded that future developments in nuclear theory and data may suggest the inclusion 
of further parametric uncertainties, most notably those related to deformation and to low-lying (3 + strengths. 

Indeed, a reliable description of the low-lying j3 + strengths is a challenging task, which calls for some improvement 
of the QRPA calculations. In fact, the 2i//3/3-decay matrix element used to fix the value of g pp is dominated by 
contributions of low-lying states of the intermediate nucleus. A recent study [19j | has shown that a better agreement 
for contributions of low-lying states to the 2^/3/?-decay matrix element can be achieved by adjusting the single-particle 
energies so as to reproduce experimental occupation numbers of neutron and proton valence orbits in 76 Ge and 76 Se. 
For a systematic analysis of this kind of effects one needs more experimental data (measuring the neutron and proton 
occupancies in particle adding and removing transfer reactions [2(j, measuring the beta strength distributions in 
charge-exchange reactions [2l|, etc.) and further theoretical studies, which go beyond the scope of this paper. 

Our results must thus be interpreted as an attempt to quantify conservatively the role of known QRPA uncertainties, 
which does not exclude that further corrections may be required by future developments in this evolving field of 
research. 

B. Comparison with other estimates and discussion 

In Fig. 2 (in the same coordinate planes of Fig. 1) we show our error ellipses at 1, 2 and 3 standard deviations 
(A% 2 = 1, 4, and 9, respectively), and superpose the latest QRPA results from Ref. [22| (dots) and the latest shell- 
model results from Ref. [23|,[24| (stars, for the available nuclei). For each nucleus, these independent rji evaluations fall 
within our estimated 3cr range, r?? ± 3(7^. Joint estimates of (rji, rjj) for couples of nuclei appear to be roughly aligned 
along (or parallel to) the major axis of each ellipse, providing an independent confirmation of positive correlations 
between the NME. The joint estimates also fall within our 3cr ellipses in most cases, with a few moderate exceptions 
in some panels of Fig. 2. We refrain, however, from enlarging our errors (or decreasing their correlations), in order to 
accommodate these few outliers within each 3cr error ellipse. A motivated revision of our estimates should be based 
on a detailed comparison of our probability distributions with analogous ones from independent calculations — rather 
than with a few sparse points from the published literature. 

Therefore, it would be useful if other theoretical groups in the Qvj3(3 field could also present "statistical samples" of 
NME calculations, as suggested in this work, so as to provide independent estimates of (co)variances for their NME 
estimates. In fact, our (co)variances cannot be directly applied to other NME evaluations which, in general, do not 
share the same set of error sources. In any case, we stress that our evaluation of QRPA uncertainties is conservative 
enough to cover the most updated, independent NME calculations within ±3tr for each individual nucleus — see 
Table VI in the Appendix. Further work is clearly needed to achieve a better convergence among the central values 
estimated in different models and, possibly, to reduce their associated errors. 

Some final remarks are in order. As already mentioned, the high correlation in each panel of Fig. 1 is mainly due to 
the fact that, if the s.r.c. model or the gA parameter are varied, all NME tend to either increase or decrease jointly. 
However, the assumption that gA is the same in all nuclei may be too strong, as the amount of quenching might 
change in different nuclei. In particular, we have shown in [25j that, by using more data besides 2vf3(3 as additional 
constraints, the fitted values of gA is not necessarily constant. Independent variations of gA in different nuclei would 
generally weaken the correlations in Fig. 1. Similarly, nucleus-dependent deformations (ignored in this work) might 
lead to a further spread of errors and to weaker correlations. In general, for any two given nuclei, the more different 
their physics (in terms of g ppi gA, deformation, etc.), the weaker their correlation (in terms of nuclear matrix elements). 
Our estimated correlations might thus be lowered in the future, should the standard assumptions in QRPA modeling 
be relaxed in different ways for different nuclei. Despite all these caveats, our p t j matrix represents at least a first, 
approximate attempt to quantify existing correlations of theoretical uncertainties. Neglecting p y - altogether would 
definitely lead to worse approximations. 
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TABLE II: 


Best current limits on 


half-lives at 90% C.L. (T< > T 90 and r< 


> t, 90 ) for different nuclei i, from [2r3|. 




% 


T 90 lv 

J i /y 




_CjXp6riII16irt 


rvei. 


76 Ge 


1.6 x 10 25 


25.204 


IGEX 


[27] 


82 Se 


2.1 x 10 23 


23.322 


NEMO-3 


[28] 


96 Zr 


8.6 x 10 21 


21.934 


NEMO-3 


M 


100 Mo 


5.8 x 10 23 


23.763 


NEMO-3 


[281 


116 Cd 


1.7 x 10 23 


23.230 


Solotvina 


[29] 


128 Te 


7.7 x 10 24 


24.886 


Geochem. 


[30] 


iso Te 


3.0 x 10 24 


24.477 


CUORICINO 




136 Xe 


4.5 x 10 23 


23.653 


DAMA 


[32] 



IV. APPLICATIONS 

In this Section we apply the previous results to cases of practical interest, in order of increasing complexity. 

A. An application not involving correlations 

As a first application (not involving correlations), we translate 90% C.L. limits on half-lives into 90% limits on the 
Majorana neutrino mass. We remind that a two-sided 90% C.L. range corresponds to ±1.64a (Ax 2 = 2.7); therefore, 
the claim in Eq. (fT4")) corresponds to 

= 25.355 ±0.118 (90% C.L., i = 76 Ge) . (15) 
and thus to the following 90% C.L. range for \i [as given by Eq. (|10[)]: 

/x ± 1.64(5 = ha- rf ± 1.64 fli ) - ± 1.64 a,) 

= -0.554 ± 0.208 (90% C.L., i = 76 Ge) , (16) 

where the two errors (st/2 and 0%) have been added in quadrature, being uncorrelated. The corresponding preferred 
range for the Majorana neutrino mass is: 

m w /eV = [0.17, 0.45] (90% C.L., i = 76 Ge) . (17) 

The best one-sided 90% C.L. limits for various nuclei have been recently reviewed in 26], in terms of half-lives at 
90% C.L. (ri > rf°), as reported in Table II. It is worth noticing that, if former data from the Heidelberg-Moscow 
experiment were interpreted as a limit on (rather than a signal of) 0^2/3 decay, the 90% C.L. bound on the 76 Ge 
half-life would be 1.9 x 10 25 y slightly stronger than the one placed by IGEX [27[ in Table II. 

The information in Table II can be transformed into 90% C.L. limits of the form [i < /i 90 via the relation 

< l^-rf^-^ + lMa^n 90 , (18) 

where we have linearly added two one-sided limits at 90%: an experimental one(— rf°/2) and a theoretical one 
( — Vi + l-64cTj). In the absence of more detailed information about the (unpublished) probability distribution of 
experimental r^'s, this is the most conservative choice. 

Figure 3 shows the results of this exercise, in terms of m^^/eV = 10 M . The shaded band on the left corresponds 
to the 90% C.L. range in Eq. (fl~5|) . while the bands on the right are obtained by inserting the rf° limits of Table II 
into Eq. (|18p . except for the very weak limit from 96 Zr which is out of scale. No experiment appears to have probed 
the 90% C.L. range preferred by the Klapdor et al. claim, although IGEX and CUORICINO have almost reached its 
lower end. 
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It is affirmed in [3l[ that the CUORICINO limit probes part of the Klapdor et al. range in mpp, seemingly in 
contrast with our results. However, the arguments in [3l| involve a comparison of two different confidence levels, 
namely, the 90% C.L. limit from 130 Te versus the 99.73% C.L. range (±3cr) from 76 Ge. The latter range is a factor 
of 3<7/1.64er = 1.83 wider than the appropriate 90% C.L. range used in Fig. 3 (left side), and thus leads to more 
optimistic conclusions. In Ref. [lj] the comparison was consistently made at the same C.L. for both nuclei, but it 
involved an intermediate step where correlations were not taken into account (see next subsection), leading again to 
an optimistic impact for the CUORICINO limit. The analysis proposed in this work shows that, actually, neither 
IGEX nor CUORICINO exclude fractions of the range claimed in [16| at comparable confidence levels, as far as our 
estimates for r\i = 77° ± Uj (and p^) hold. 



B. Comparison of half-lives in a couple of nuclei 



Here we consider a more direct comparison via observable half-lives in two nuclei, bypassing the unobservable 
Majorana mass mgg. We take two different nuclei i and j, characterized by nuclear matrix elements 77^ = 77? ± Oi and 
r\j = r) j ± Oj with correlation pij. A positive Ov/3/3 signal in the first nucleus (7^ = t® ± Sj) translates into a favored 
range for the second nucleus (73 = rj* ± Sj) as follows. 

From Eq. (jTOJ) one obtains, by difference, 



where 



Aij ± eij , (19) 



A 4 , = r? - r° = ( 7j - 7 + 2(*?° - „?) , (20) 



the error €ij being obtained by summing in quadrature the correlated uncertainties associated to the difference 

2(1?-^ 

4 = 4 ( ct ? + °j - tpwi) ■ (21) 



Note that, if the correlation term —2pij(7i<7j were neglected, Cy would be overestimated. The error Sj associated to 
Tj = rf + is obtained by summing in quadrature the uncorrelated errors Sj and e,j , 



a? = a) + 4 . (22) 
As a result, the error sj has a nonzero correlation with the error Sj, as given by Sj Sj — s|, namely, 

r« = ^ . (23) 

If we apply the above results to i = 76 Ge and j = 130 Te, then the claim by Klapdor et al. in Eq. [T?l n = 
25.355 ± 0.072, implies that Tj = 24.786 ± 0.161, with error correlation nj — 0.447. Figure 4 shows the corresponding 
error ellipse at 1.64<r (90% C.L.), in the plane charted by the Qvf}(3 half-lives of the nuclei i = 76 Ge and j = 130 Te. 
The ellipse can be thought as the combined result of two independent constraints, shown as 90% C.L. bands. The 
horizontal band corresponds to the experimental claim Ti — 25.355 ± (1.64 x 0.072). The slanted band corresponds 
to the theoretical limits placed by our QRPA estimates on the ratio Tj/Ti, namely, Tj — Ti = A ± (1.64 x ey). 
Note that the projection of the ellipse on the x-axis provides the range preferred at 90% C.L. for the 130 Te half-life: 
Tj/y — [0.33, 1.12] x 10 25 . Projections for other nuclei can be similarly derived, as reported in Fig. 5. 

Figure 5 shows the two-sided ranges preferred by the Klapdor et al. claim at 90% C.L. (shaded rectangles on the 
right), as well as the one-sided 90% C.L. limits from Table II (bands on the left), for the same nuclei as in Fig. 3. 
The two-sided limits involve the use of NME errors and correlations, except for 76 Ge, which is a purely experimental 
input. Once more, we see that none of the existing limits can exclude a fraction of the range favored by Klapdor et al. 
fig at a comparable confidence level, alth oug h IGEX and CUORICINO have almost reached it. The more optimistic 
claim about the CUORICINO impact in [31( was based on a larger favored range for the 130 Te half-life, as obtained 
by ignoring correlations in the estimate of Eq. (j!?lj) . 

We emphasize that the contents of Figs. 3 and 5, although similar, are not equivalent. The comparison of ex- 
perimental sensitivities in Fig. 3 is made in terms of a derived quantity (rripp), while in Fig. 5 it is directly made 
in terms of observables (Ti). One-sided bounds in Fig. 3 are obtained by linearly adding 90% C.L. theoretical and 
experimental limits [Eq. I|18[)]. while in Fig. 5 only the latter limits are used; conversely, theoretical errors are used 
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TABLE III: Prospective half-life sensitivities at 90% C.L. (Xf°) for different nuclei i in promising future projects, as reported 
in 0. 

i Tf°/y Project 

76 Ge 2.0 x 10 26 GERDA, MAJORANA 

82 Se 2.0 x 10 26 SuperNEMO 

130 Te 2.1 x 10 26 CUORE 

136 Xe 6.4 x 10 25 EXO 



with full correlation information in the allowed (two-sided) bars of Fig. 5. We think that a comparison in terms of 
observables, as in Fig. 5, provides a more faithful representation of the current §v(3(3 decay sensitivities. 

We conclude this subsection by discussing the 90% C.L. prospective sensitivities (in terms of TA of the most 
promising future 0^/3/3 projects. Table III reports such limits, according to the recent review in Ref. [26|]. The values 
in Table III are largely beyond the two-sided favored ranges in Fig. 5, except perhaps for 136 Xe, where the expected 
sensitivity is only a factor < 2 beyond the Klapdor et al. favored range. This gain may be insufficient if one requires a 
more demanding check of the claim, at a confidence level significantly higher than 90%. It should be added, however, 
that all the projects in Table III expect to proceed in a second phase of operation with larg er exposures and lower 
backgrounds, improving the quoted sensitivities by, possibly, another order of magnitude [26j. 



C. Combination of half-life data from several nuclei, and degeneracy effects 

Let us consider a future, optimistic situation where §v(3(3 decay is established in N different nuclei, with measured 
half-lives 

Ti = T°± Si (i = l,...,N) . (24) 

Assuming that Qv[3(3 decays proceed only through light Majorana neutrino exchange, these measurements will fix one 
unknown parameter (/i) via a set of N linear equations analogous to Eq. (p~0|) . 

r° ± Sl = 7l - 2(?7? ± (Ti) - 2fi (i = 1, . . . , N) , (25) 

where the experimental errors Sj are, in general, uncorrelated (being obtained in independent experiments), while the 
theoretical errors Ci have nontrivial correlations pij (being obtained within the same QRPA model). 

This overconstrained system can be solved by the least-squares method, i.e., by minimizing the \ 2 function 

X 2 (M) = E( T ? - ^ + + 2 v) w v (rf ~ 7i + 2*7° + 2 M ) , (26) 

where the weight matrix Wij is the inverse of the total covariance matrix (including experimental and theoretical 
errors), 



\W)ij = SijSiSj + ApijCTicrj . (27) 



The x 2 function is quadratic in p, 



where 



X 2 (m) =ap 2 +bp + c , (28) 

4£Wy, (29) 

ij 

4^W ij (rP- 7i + 2 ?7 ?) , (30) 

ij 

- 7l + 2^)^(7-? - 7 , + 2r?°) . (31) 
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TABLE IV: Combination of any among the four hypothetical half-life data Tj in Eq. (|35f) with experimental uncertainty 
STi/Ti = 20%. Results are given in terms of the total lcr error S on the parameter n = log 10 (m/3/3/eV), including theoretical 
uncertainties without and with correlations. Bullets indicate the data included in the evaluation (from 1 to 4 data). 



# of data 


76 Ge 


82 Se 


130 Te 


136 Xe 


5 (w/o corr.) 


S (with corr.) 


1 











0.128 


0.128 


1 










0.141 


0.141 


1 










0.163 


0.163 


1 








• 


0.191 


0.191 


2 











0.095 


0.128 


2 











0.100 


0.128 


2 










0.106 


0.127 


2 




• 


• 




0.107 


0.141 


2 




• 




• 


0.114 


0.141 


2 






• 


• 


0.124 


0.163 


3 


• 


• 


• 




0.082 


0.127 


3 


• 


• 




• 


0.085 


0.127 


3 


• 




• 


• 


0.089 


0.127 


3 




• 


• 


• 


0.093 


0.140 


4 


• 


• 


• 


• 


0.075 


0.127 



The minimum value Xmin an< ^ the one-sigma shift Xmin + 1 are reached for \x — \x§ and /i = [1q ± 5, respectively, where 

Mo = -I , (32) 
* = 4= . (33) 

Xmin = C ~ ^ • (34) 

The fit is acceptable if X m t n /(iV — 1) ~ 1. Much higher value of Xmin might signal, e.g., new physics beyond the 
standard mechanism of §v(3(3 decay via light Majorana neutrinos (barring experimental and theoretical mistakes). 
However, the analysis of nonstandard mechanisms is beyond the scope of this work. 

As a practical example for the standard §v(3(3 case, we consider decay searches in each of the four nuclei reported 
in Table III, in the hypothesis that the true value of mpp is 0.2 eV (i.e., fj, = —0.70), close to the lower end of the 
range in Eq. (fT7|) . We assume that the experiments will measure the expected values for the half- lives Tj, 

4.43 x 10 25 ( 76 Ge) , 

m , ,/VY = 0.2 = TVv = { 1M X W Z 5!L Se) ' (35) 

' 1.20 x 10 25 ( 130 Te) , 

3.43 x 10 25 ( 136 Xe) , 

with a fractional uncertainty STi/Ti — 20% (corresponding to Sj = 0.08). By construction, the best fit to any 
combination of these mock data gives back /iq = —0.7 and X m in~0- The relevant output parameter is then the 
reconstructed fi uncertainty, S, from Eq. (f33|) . 

Table IV shows the 8 values, for all possible combinations of mock data from the four nuclei (ranging from a single 
nucleus to all of them). We comment first the results in the 6th column, which are obtained by (incorrectly) switching 
off correlations, i.e., by setting pij = 8ij, as it is often done in the literature. 

Without correlations, the error S is given by the familiar combination of total errors from independent data, 

^ = ?a 2 + ( Sl /2) 2 "^^' (36) 



where we have used the fact that any of the <7j is a factor of 3-4 greater than s,/2 = 0.04. Although the error 6 is 
dominated by theoretical uncertainties, it decreases by increasing the data sample (see 6th column of Table IV), as 
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a consequence of (incorrectly) assuming no correlations. Formally, the combination of all the four data would then 
provide the estimate fi — — 0.7 ± 0.075, corresponding to nipp ~ 0.2 ± 0.035. 

Unfortunately, including correlations spoils this nice result. Table IV (last column) shows that, with good ap- 
proximation, the uncertainty S cannot be much better than the smallest theoretical uncertainty Oi among the set of 
nuclei included in the fit. Indeed, even with all four nuclei one obtains S — 0.127, nearly the same as 5 = 0.128 from 
the single nucleus 76 Ge (characterized by the smallest theoretical error, a = 0.122). Therefore, regardless of how 
many accurate experiments are combined, the final accuracy for our test-case Majorana mass will not be better than 
// ~ -0.7 ± 0.13, namely, mpp ~ 0.2 ± 0.06. 

The degeneracy effect induced by correlations can be easily understood in the limiting case of equal and completely 
correlated theoretical errors (e^ = a and pij = Sij). In this case, the QRPA uncertainties would reduce to a common 
shift rji — » rji + S for all nuclei, where 5 6 [— cr, +<j] within one standard deviation. [In Fig. 1, the ellipses would 
collapse to "segments" with 45° slope in all panels.] A common shift of all rji is degenerate with a shift /_t — > fi — 6 via 
Eq. (Hi. 

n = 7i - 2(77, + S) - 2(ji - 6) , (37) 

and, thus, the parameter fx is affected by an irreducible uncertainty S — a. For unequal NME errors ov, the most 
accurate one dominates in equations like Eq. (f3"T)) and thus 

S ~ min{(Tj} , (38) 

as anticipated. 

The difference between S estimates without or with correlations, in the combination of data from N different 
nuclei, is striking. Without correlations, and for comparable theoretical uncertainties, the error S would scale as yN 
[Eq. (f3"6"| ]. Including correlations, the error 5 becomes dominated by the single, most accurate NME, irrespective of 
N [Eq. (|38|) ]. One should thus reduce not only the size, but also the correlations of theoretical errors, in order to fully 
exploit the mpp sensitivity of future, multiple-isotope Ov/3/3 searches. 



D. Prospective constraints on the absolute neutrino mass and Majorana phase 

The Majorana mass in Of/3/3 decay {mpp) is one of the most sensitive probes of the absolute neutrino mass scale m„, 
together with the effective neutrino mass in beta decay (mpp) and the sum of the three neutrino masses in cosmology 
(S); see [l4[ for updated bounds. It is tempting to combine prospective data on (mpp, mp, S) in the optimistic case 
of a possible sig nal "waiting around the corner", i.e., for masses close to the current conservative cosmological bound 
E < 0.6 eV Hi SI HI: 

mi ~ ?7i2 — T7.3 = m-v — 0.2 eV . (39) 



For the sake of simplicity, within current neutrino oscillation phenomenology 1J| , we approximate the mixing matrix 
values as: 

U 2 el ~ 0.69 , (40) 

t/ e 2 2 ~ 0.31 e** , (41) 

U 2 e3 ~ , (42) 

where cj> is an unknown Majorana phase. For nearly degenerate masses it is thus [l2| 

m ~ m v , (43) 

S ~ im v , (44) 

mpp ~ m u f , (45) 

with 

/- \U 2 el +U 2 e2 \ e [0.38, 1] , (46) 

where the upper (lower) end of the range is obtained for the CP-conserving case e 1 ^ — +1 (e 1 ^ = —1). 

Let us test the above scenario with mock data, having the following central values and fractional la errors: 

to^ ~ 0.2(1 ±0.5) eV , (47) 

E ~ 0.6(1 ±0.3) eV , (48) 

m Pf ) ~ 0.2(1 ± 0.3) eV . (49) 
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In the above equations, the 50% uncertainty on mp corresponds to the smallest la error estimated for the upcoming (3- 
decay experiment KATRIN {Smp ~ 0.1 eV) [36| . A 30% uncertainty on E seems appropriate (and even conservative) 
for a signal in next-generation cosmological data [53, [38|. The putative 30% uncertainty on mpp reflects the discussion 
in the previous subsection. 

Combining the "data" in Eqs. (147)) and (|48|) . one obtains 

m v ~ 0.2(1 ± 0.25) , (50) 

which, together with Eq. (j45|) and the "datum" in Eq. (|49| . imply 

/~1±0.4. (51) 

This result, compared with the range in Eq. (|46|) . would slightly prefer one CP-conserving case (e 1 ^ = +1) over the 
other {e 1 ^ = —1), at the level of ~1.5er. Therefore, in an optimistic — but not completely unrealistic — scenario with 
degenerate neutrino masses, such as the one considered above, a possible determination of m„ ~ 0.2 eV with ~25% 
accuracy (via 771,3 plus E) might be accompanied by some indications about the Majorana phase <f> (via Tripp). In this 
sense, we feel sympathetic towards more encouraging viewpoints [3!|, H(| than those expressed by a "no-go detection" 
for <p [4lj |. although a real "measurement" of tf> remains undoubtedly very challenging, even in the most favorable 
scenarios. 



V. SUMMARY AND PROSPECTS 

Nuclear matrix elements for Qvj3(3 decay are affected by relatively large theoretical uncertainties. Within the QRPA 
approach, we have shown that, within a given set of nuclei, the correlations among NME errors are as important as their 
size. We have made a first attempt to quantify the covariance matrix of the NME, and to understand its effects in the 
comparison of current and prospective Qvf}f3 results for two or more nuclei. The effects have been clarified through a 
series of examples, involving an increasing number of observables. It turns out that correlations may severely limit the 
accuracy in the reconstruction of mpp from any number of 0v(3f3 observations in different nuclei, due to a degeneracy 
between NME and mpp uncertainties. In particular, the fractional error on mpp is ultimately dominated by a single 
fractional NME uncertainty (the smallest one, among the set of nuclei considered). Breaking correlations between 
different nuclei is thus an important goal, which requires constraining (and improving) the theoretical model of each 
nucleus by means of many independent data (not only 2v(3[3 data as currently used). In this way, systematic effects 
common to all nuclei may be reduced. Another relevant goal is to compare correlation estimates in future independent 
calculations (e.g., QRPA versus shell- model) . While pursuing such a long-term theoretical and experimental program, 
a covariance analysis like the one proposed in this work may represent a useful tool, in order to correctly estimate 
current or prospective sensitivities to 0v(3f3 decay and to Majorana neutrino parameters. 
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TABLE V: Comparison of Gim 2 estimates (in units of 10 15 y _1 ) for qa = 1.25. The second column refers to the calculations 
reported in [43t ] for ro = 1.1 fm, as used in this work. The third and fourth columns refer to independent estimates [il. EH. I46T] 
for ro = 1.2 fm, rescaled by a compensating factor /q = 1.2. 



Nucleus Ref. [43J Refs. [44, 45] Ref. [46] 

76 Ge 7.93 7.67 7.57 

82 Se 35.2 33.8 32.8 

96 Zr 73.6 70.2 68.4 

100 Mo 57.3 54.8 52.8 

116 Cd 62.3 59.3 56.2 

128 Te 2.21 2.20 1.99 

130 Te 55.4 53.2 49.7 

136 Xe 59.1 56.8 52.4 



TABLE VI: Estimates of r\i — log 10 \M-\ for each nucleus, as derived from the recent QRPA calculations in 123] (see Tab. 1 
therein) and shell-model calculations in [24J (see Tab. 7 therein) after appropriate rescaling, in order to match the conventions 
used in this work. The estimates of [24| refer only to a subset of nuclei and to qa = 1.25. The s.r.c. used (Jastrow or UCOM) 
are explicitly reported. The last two rows report the upper and lower ends of our three-standard-deviation ranges r/® ± 3<Ti (for 
any s.r.c. and qa), which embrace all the above r\i estimates. 



Ref. 


s.r.c. 


9A 


76 Ge 


82 Se 


96 Zr 


100 Mo 


116 Cd 


128 Te 


130 Te 


136 Xe 


[22] 


Jastrow 


1.00 


0.471 


0.313 


0.261 


0.312 


0.331 


0.396 


0.374 


0.222 


[22] 


UCOM 


1.00 


0.582 


0.427 


0.400 


0.451 


0.435 


0.531 


0.501 


0.335 


[22] 


Jastrow 


1.25 


0.564 


0.401 


0.274 


0.396 


0.441 


0.488 


0.435 


0.271 


[22] 


UCOM 


1.25 


0.687 


0.529 


0.452 


0.553 


0.554 


0.639 


0.584 


0.406 


[24] 


Jastrow 


1.25 


0.320 


0.297 








0.328 


0.285 


0.204 


[24] 


UCOM 


1.25 


0.407 


0.380 








0.418 


0.382 


0.299 


This work 
This work 


Lower limit at 3<r level 
Upper limit at 3a level 


0.269 
1.001 


0.166 
0.976 


-0.703 
0.779 


0.017 
0.989 


-0.046 
0.854 


0.072 
0.996 


0.024 
0.972 


-0.307 
0.815 



APPENDIX 

This Appendix clarifies the role of different conventions about the 0vj3j3 phase space factor Gi, the axial vector 
coupling §a, the nuclear matrix elements M' tl and the nuclear radius i? nU ci, in comparison with other authors. An 
agreement on common conventions would be desirable in the future, to avoid possible confusion or ambiguity (see 
also [13). 

According to usual definitions, the phase space Gi contains a factor (g 2 A / Rnuci) , where i? n uci = tqA 1 ^ . In this work, 
the adopted values of Gi refer to ro = 1.1 fm and gA — 1-25 [43fl. while changes of gA are conventionally embedded in 
M[ (rather than in Gi) via the prefacto r (ci a/1.25) 2 in Eq. ^ pLIlljj]. In order to match such convention, alternative 
calculations of Gi using ro = 1.2 fm [iilEaE^] must be rescaled by a factor /g ~ 1.2 (where /o — 1.1 — 1.2/1.1) [42| . 

Table V compares three different phase-space calculations (in terms of Gi ra 2 ,), all normalized to the same reference 
values gA = 1-25 and r = 1.1 fm. One can notice residual differences of ~ 5% between the results of [43[ and [11, El], 
and of ~ 10% between those of [43| and [46[ , presumably due to different approximations used to evaluate the electron 
wave function and the screening corrections. In our opinion, a typical uncertainty for the computed Gi values may 
be estimated as ±5%, corresponding to a variation Sji ~ ±0.02 for the 7, values in Table I. Such minor error, being 
much smaller than the theoretical and experimental uncertainties considered in this work (07^ <C &i, Sj), has been 
ignored — but it might become more important in the future. 

Concerning the nuclear matrix elements |M||, the values calculated in [22j (QRPA) and 0, [24[ (shell model) refer 
to ro = 1.2 fm, and must be rescaled by a factor l//o for comparison with the NME used in this work. Furthermore, 
since the values in [22j do not embed the prefactor (<ju/1.25) 2 , they must be rescaled by another factor 1.25~ 2 in 
the subcase gA = 1; this further rescaling is not necessary for the NME values in [24|. Table VI reports the rescaled 
values of \M- \ from [22j and [24| (in terms of logarithms r/i), as also used in Fig. 2. These r/i values are all contained 
within our estimated three-standard-deviation ranges 77° ± 3<Ti, which are reported in the last two rows of Table VI. 
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FIG. 1: Scatter plot of estimated QRPA values for the (logarithms of) nuclear matrix elements (r/i, rjj) for each couple of nuclei 
(i, j), together with the error bars induced by g pp uncertainties. In each panel, also shown is the la error ellipse, conservatively 
estimated on the basis of the scatter plots. See the text for details. Color code for s.r.c: blue (Jastrow), red (UCOM). 
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FIG. 3: Range of allowed at 90% C.L. by the 0vf3f3 claim of [l6| ]. compared with the 90% limits placed by other 

experiments. The comparison involves the NME and their errors, as estimated in this work. 
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FIG. 4: Example of theoretical and experimental constraints at 90% C.L., in the plane charted by the Of/3/3 half- lives of 
76 Ge and 130 Te. Horizontal band: range preferred by the Ovf3f3 claim of [l6l ]. Slanted band: constraint placed by our QRPA 
estimates. The combination of the two constraints provides the shaded ellipse, whose projection on the abscissa gives the range 
preferred at 90% C.L. for the 130 Te half life. 



18 



Exp. limits, 90 % C.L H Klapdor et al. ( + NME), 90 % C.L 




FIG. 5: Range of half-lives Ti preferred at 90% C.L. by the 0vf3f3 claim of fl6j | . compared with the 90% limits placed by other 
experiments. The comparison involves the NME and their errors, as well as their correlations, estimated in this work. 



